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ABSTRACT 

Long-range, over-the-horizon (transhorizon) radio wave propagation is considered 
for the case of the moon. In the event that relay satellites are not available or otherwise 
unwarranted for use, transhorizon communication provides for a contingency or backup 
option for non line-of-sight lunar surface exploration scenarios. Two potential low- 
frequency propagation mechanisms characteristic of the lunar landscape are the lunar 
regolith and the photoelectron induced plasma exosphere enveloping the moon. 

Although it was hoped that the regolith would provide for a spherical waveguide which 
could support a trapped surface wave phenomena, it is found that, in most cases, the 
regolith is deleterious to long range radio wave propagation. However, the presence of 
the plasma of the lunar exosphere supports wave propagation and, in fact, surpasses the 
attenuation of the regolith. Given the models of the regolith and exosphere adopted here, 
it is recommended that a frequency of 1 MHz be considered for low rate data 
transmission along the lunar surface. It is also recommended that further research be 
done to capture the descriptive physics of the regolith and the exospheric plasma so that a 
more complete model can be obtained. 

This comprehensive theoretical study is based entirely on first principles and the 
mathematical techniques needed are developed as required; it is self-contained and 
should not require the use of outside resources for its understanding. 
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I. Introduction 


The subject of this memorandum is a theoretical study of the possibility of long 
range electromagnetic wave propagation beyond the horizon (i.e., transhorizon) for 
purposes of contingency communications to support lunar surface exploration activities 
when the use of a relay satellite is not possible or warranted. For a source of radiation on 
or above the lunar surface, diffraction will of course be a mechanism for such 
transhorizon propagation as it is on the Earth. This effect is prevalent only at low radio 
frequencies in the KHz or low MHz range. Thus, lunar transhorizon propagation can 
only be considered at low frequencies and hence, low data rates. However, the lunar 
environment has the features that can also support long-range propagation. 

The moon possesses two characteristic features that can potentially lend 
themselves to the support of transhorizon radio frequency propagation along its surface. 
First, the moon has a regolith [1], i.e., a thin fragmental layer of rock and minerals as well 
as glass, that covers all the rock fonnations on the moon. Its thickness subtends 10 
meters to 100 meters below the surface; the thickness most likely varies along the surface 
especially near crater rims and mountain ranges. The source of the regolith is meteorite 
bombardment and related secondary processes as well as volcanic eruptions. 

If the electrical characteristics of the regolith are sufficiently different than those 
of the core, the regolith has the potential to act as a spherical waveguide for low 
frequency electromagnetic radiation within which trapped surface waves or “whispering 
gallery” waves can travel to great distances along the surface. 

In addition to the regolith, the moon also possesses an “ionosphere”, or more 
appropriately, an exosphere, a layer of electrons that is apparently formed by the 
liberation of electrons from the lunar surface by the solar wind. This photoelectron layer 
has been observed by the Apollo Lunar Science Experiment Package (ALSEP) [2] as 
well as by the Russian Luna 22 spacecraft [3], It is important to note that this 
photoelectron layer is not a thermal ionospheric plasma like that of the Earth, which 
results from photoionozation of neutral exospheric molecules [4]. As a result of the 
photoelectron source of the lunar exosphere, this phenomenon can only be expected to 
occur on the dayside of the moon 

As will be displayed below, the charge density of the lunar exosphere begins at 
the surface and increases up to the height of 10 km after which it decreases again. In 
such a scenario, electromagnetic waves, with a source on the lunar surface, will not 
reflect beyond the horizon as in the case of the layered ionosphere of the Earth, but will 
again be guided around the surface of the moon as in the case of a dielectric waveguide 
of varying refractive index. 

The purpose of this work is to theoretically study the effect of the lunar regolith as 
well as its exosphere plasma on low-frequency radio wave propagation. Thus, there are 
two propagation mechanisms that will be collectively studied to determine the potential 
of long-range, transhorizon radio wave transmission, i.e., that of the regolith and the 
contribution of the exosphere. Section II will detail the electrical properties of the 
regolith and the exosphere. In both cases, the properties have only been so far remotely 
obtained and can, over the course of future lunar exploration, change. The mathematical 
model of these propagation mechanisms will be kept general enough to accommodate the 
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evolution of the quantitative understanding of the regolith and exosphere. Section III will 
present the mathematical model of the lunar propagation environment. Propagation 
equations will be established for the three regions of the core, regolith, and the exosphere. 
These differential equations will then be solved by matching boundary conditions across 
the regions. Two cases will then be considered, i.e., the case where there is only a 
regoilth with no exosphere and the case in which the exosphere is present. This exercise 
will allow for the development of the mathematical technique using the simpler case of 
just the regolith and, in addition, will allow one to isolate the effect of the regolith alone 
on transhorizon propagation. As far as mathematical technique is concerned, every effort 
is made to make this study self-contained and mathematical detail will not be spared. 

Two appendices are also provided to help with the rigor. 


II. THE CORE, THE REGOLITH, AND THE EXOSPHERE 

Figure 1 displays the model of the lunar propagation environment as envisioned in 
this work. Here, the Moon is modeled to possess a core (Region 1) of dielectric 
permittivity s x followed by a very thin regolith (Region 2) of thickness 8 and permittivity 
s 2 \ both the core and regolith give rise to an overall lunar radius a. 



Figure 1. Defining the Propagation Regions of the Moon 
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Here, g 2 and g, are ‘local’ radial coordinates measured, respectively, from the surface of 
the core and the surface of the regolith. Later, a spherical coordinate system will be 
introduced in which the radial coordinate of which will be reckoned to the center of the 
moon. 

The outer surface of the Moon is taken to have an exosphere (Region 3). The electron 
concentration profile within this region gives rise to a tenuous plasma of permittivity 
profile e 3 (£ 3 ). 


a. Permittivity and Conductivity of the Core and Regolith 

A majority of the electrical characteristics of the core and regolith are, of course, 
indirectly obtained [5,6] and their estimates seem to vary. In a related investigation [7], 
the moon is taken to have a regolith that subtends the first 100 meters below the surface, 
after which the core is reached. However, as mentioned above, some studies indicate that 
the regolith thickness can be a minimum of 10 m. In this study, two values for the 
thickness of the regolith will be considered, viz., 8 = 50 m and 8 = 100 m. The 
pennittivity of the core is taken [7] to be within the range from s x = 5 to s x = 1 0 with an 

associated conductivity within the range from <j x = 10 4 mhos/m to a x = 1 0 mhos/m. 
Similarly, the regolith region is characterized [7] by a pennittivity range from s 2 = 2 to 
£, = 3 and a conductivity range from cr 2 = 10 mhos/m to a 2 = 10 mhos/m. For the 
purposes of the present investigation, the values of the relevant parameters are taken as 
displayed in Table 1. 


REGION 

PERMITTIVITY 

CONDUCTIVITY 

Core 

s x = 5 

<r x = 5 * 10~ 4 mhos/m 

Regolith (8 =50 m, 
8 = 100 m) 

^ 2 =3 

<j 2 = 5 * 1 0 mhos/m 


Table 1. Electrical Parameters of the Moon 


b. Permittivity and Conductivity of the Exosphere 

As mentioned above, the lunar exosphere has been experimentally observed [3], 
the properties of which are shown in Figure 2 which displays the height profile of the 
electron concentration of the lunar exosphere. As can be seen form the figure, the 
electron concentration has a maximum N m = 1000 cm 3 occurring at the height 
^ 30 = 10 km. Here, is the ‘local’ coordinate measured from the lunar surface, as 
shown in Fig. 1 . The electron concentration profile of the exosphere is most easily 
modeled as a parabolic profile with its apex at g x = g V) . 
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Figure 2. Electron Concentration Profile of Lunar Exosphere 


Thus, letting rj be a constant, one can write for this concentration 

^3) = A^-^3-^3o) 2 C 1 ) 

This constant is determined from the fact that the concentration at the surface N(0) is 
given, by Eq.(l) as 


jv (°) = 'V,,-')(‘5 jo) 2 (2) 

Hence, I] = AN /c,l } where AN = N m - N(0) . From Fig.2, the differential concentration of 
the electrons is AN = 400 cm 3 . Thus, the model for the electron concentration profile of 
the lunar ionosphere is given by 




f 

-AN ■ 
V 


lizS 

^30 


30 


v 

) 


(3) 


The associated permittivity profile of this ‘p ara t>olic’ exosphere is found by 
applying a well-known model for free electron plasma behavior in an electric field [8] 
with e~ ,m time dependence that gives 


e 3 (£,) = ! 


Ane 2 N (£ 3 ) 
mO)[o) + ) 


(4) 


where e and m are, respectively, the charge and mass of an electron, (0 is the angular 
frequency of the electromagnetic radiation and v eff is the effective collision frequency for 

the electrons. This equation can be re-written as 
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^fe) =1 


( 5 ) 


Ane 1 N(^ i ) An 

( 2 \ \ +? — a : 
m[a> +\r ff ) co 



where 


^3 



m (oj 2 + vl ff ) 


( 6 ) 


is the conductivity of the exosphere plasma. However, in a tenuous plasma as is 
represented by Fig.2, one has that v eff « 0 thus rendering the plasma conductivity zero 

and allowing Eq.(5), after substitution of Eq.(3) to be written as 



+ £■ 


3A 


6- 


^30 


&0 


(V) 


where 


_ 4 ne 2 N m _ Ane 2 AN 

^30 ^ 2 ’ ^3 A — 2 v®j 

mco mco 

Equations (7) and (8) constitute the permittivity model of the lunar exosphere. 

Using the values of the parameters quoted above in the text as well as 
m = 9. 1 x 10 28 g and e = 4.8 x 1CT 10 esu, Eq.(7) for £, = 0 is plotted versus the operating 
frequency as is shown in Fig. 3. The graph clearly shows that at frequencies below 220 
KHz, electromagnetic radiation does not propagate at the surface; this is the characteristic 
plasma frequency for the plasma at the surface. Only frequencies above this threshold 
can be transmitted along the lunar surface characterized by the profile of Fig. 2. Figure 4 
displays the behavior of the exospheric permittivity as a function of height above the 
surface at a frequency of 500 KHz. A minimum exists at the characteristic height . 
This interesting circumstance suggests that a ‘trough’ exists in the permittivity 
distribution along which electromagnetic waves may be guided in to a transhorizon 
location. 

Now that the regions of propagation have been identified and their prevailing 
parameter values been established, one can now begin to model the electromagnetic 
propagation behavior of these regions and how they all couple to contribute to long range 
transhorizon transmission. 
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Figure 3. Permittivity of Exosphere vs. Operating Frequency 



& ( cm ) 


Figure 4. Permittivity of Exosphere vs. Height at 500 KHz 


III. MODELING OF THE LUNAR PROPAGATION ENVIRONMENT 
a. Derivation of the Propagation Equations for the Three Regions 

As introduced above, the lunar propagation scenario is as depicted in Figure 1 . 
There are three electrically disparate regions defined by the core (region 1), the regolith 
(region 2), and the exosphere (region 3). The spherical symmetry of the situation 
naturally suggests that one consider a spherical coordinate ( r,6,(j ) ) system be attached to 
the moon, the center of which is placed at the center of the moon. Here, r is the radial 
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distance from the center, # is the colatitude, as reckoned from a polar axis which is taken 
to be along the source of the electromagnetic radiation placed on the surface, and <j) is the 
azimuthal angle. The source is taken to be a vertical electric dipole placed at a distance 
r s from the surface within region 3; such a source radiates transverse magnetic (TM) 
waves. Regions 2 and 3 are also characterized by the distances above their respective 
surfaces. Taking the radius of the moon to be a and the thickness of the regolith to be 8, 
region 2 is defined by the radial coordinate where a - 8< g 2 < a . Similarly, region 
three is defined by the radial coordinate g, where a < g 3 . 

The modeling approach must necessarily be based on the Maxwell equations that 
govern the behavior and the dynamics of the electromagnetic field. So as not to deter 
from the presentation of the subject matter, the Maxwell equations as well as their 
solution in terms of a scalar potential for TM electromagnetic waves is given in Appendix 
A. There, it is shown that the field components of such a wave, with fields that are 
independent of the azimuthal coordinate (f), are given by 


E= — 


1 d 

rsin# 36 \ 


sin # 


36 


(9) 


for the radial component of the electric field, 


E e (r,0) 


1 3( ; 3U e (r,df 
r 3r\ 36 , 


( 10 ) 


for the colatitude field component and 


H t (r,&)= +ik 0 n 2 


^A(a#) 

36 


( 11 ) 


for the azimuthal component of the magnetic field, where n 2 = s is the refractive index of 
the region in question. Here, U e (r,6 ) is the scalar potential for TM waves given by the 
solution of 


-4j( ru .) + 

r or 


1 


'd f 


V r sin#/ 


36 


sin# 


3U^ 

36) 


+ 


1 




V r sin" #A 


+ k 2 U= 0 


( 12 ) 


J 


Again, so as not to deter from the present discussion, the method of solution of Eq.(12) is 
presented in Appendix B in addition to its preferred form. Thus, the solution to Eq.(12) 
is taken to be of the form 


U. = U,(r,e)=A 


U6 


P, (- cos #) 


(13) 
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where P t (-cos 6) is the I th order Legendre polynomial in the argument -cos 6 (See 
Appendix B for an explanation of this choice rather than using the traditional +cos6f) 
and R, (/•) is given by the solution of 


d%(r) 

dr 2 



/(/+l)' 

-J- 


M>-)= 0 


(14) 


where k = ~Jak 0 is the wavenumber in the region defined by the permittivity a and 
k 0 = 2nj A is the free-space wavenumber for electromagnetic radiation of wavelength A. 

The function P } (-cos (9) is only a function of the colatitude and is therefore 
independent of which of the three regions is being considered. However, the radial 
component R,(r ) of the scalar potential is indeed a function of the region through the 
wavenumber parameter k. Thus, there are three separate versions of Eq.(14) that must be 
considered, viz., 


d 2 R h (r) 
dr 2 


+ 


1(1 + 1 ) 


2 

r J 


*,,(,')= 0, 


k(=£\kl, 0 <r<a-8 


for region 1 (the core), 


d%(r) 

dr 2 


+ 



1 ( 1 + 1 ) 


r 2 J 


R l2 (r) = 0, 


k\ = s 2 kl, a-8<r < a 


(15) 


(16) 


for region 2 (the regolith), and 


d 2 R h ( r ) 

dr 2 


+ 



1 ( 1 + 1 ) 


r 2 J 


R h (r) = 0, 


k; (g, ) = 6\ (g, )k 2 , a < r, g, = r- a 


(17) 


for region 3 (the exosphere). 


b. Solution for Region 1 

One must first consider the field distribution within the core region; this will then 
be used to calculate the impedance boundary condition at r = a - 5 to which the fields 
within the regolith can be connected. To this end, consider the solution [9] of Eq.(15) for 
the scalar potential within the lunar core 

R li (r)=r[AJ l (k 1 r)+B 1 n,(k l r)] r<a-S (18) 

where /) (k t r) and n,(k/) are spherical Bessel functions given by 
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in terms of the Bessel function of the first kind J +1/ (x) and the Bessel function of the 
second kind A |+1/ (x) (Neumann function). Hence, Eq.(13) becomes, for the core 
region, 

u e x (u °) = l A Ji ( k /)+ B i n ! ( k i r j} p i (- cos 6) (20) 

However, due to the properties of the functions j,(x) and n t (x), for the solution U (r, 0) 
to be finite at r = 0, one must have B i = 0 , 

U e x ( r > 0 ) = A Jl ( k /) P ! (- C0S 0 ) ( 2 1 ) 


The ‘matching’ of the solutions within region 1 and those of region 2 which is to 
follow, will be done through the use of the surface impedance at r = a- S. This 
impedance is defined by 


Z,= 


EjrM) 


where a R =a-S. From Eqs.(10) and (1 1), 


E «X r - 0 >-^ 

r a - 


cu.X'-.oy 


do 


h a (r,0)= +i k 0 


dU ex {r-,e) 


do 


n x = e x 


( 22 ) 


(23) 


(24) 


The use of the impedance concept of Eq.(22) is tantamount to employing the continuity 
boundary conditions of Eqs.(A48) and (A50) of Appendix A. Thus, substituting Eq.(21) 
into Eqs.(23) and (24) gives 


r 


ii ( k i r ) 


+ r 


dji ( k s) 


oP ! (-cosO) 


36 


do 


11 <h (p0)= +ik () n;A ] j l (k i r f P ' ^ ^ 


(25) 

(26) 


Using these expressions in Eq.(22) yields 
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2, 



1 | 1 dj,(k x r) 

a R ji(K r ) dr 


( i \ 1 djXk^A 

« +1 1 7 . v ’ (27) 

U 0 ejjfar) dr ^ 

The evaluation of j ,(k\d) al r = a R » I requires asymptotic expansions of this function 
which can become rather unwieldy. However, one can adopt the method of [10] and 
consider the solution to the asymptotic form of the original equation, Eq.(15) with r = a R . 
What one does in essence is to assume a thin shell about the radius a R . To this end, 
Eq.(15) is approximately written 


d 2 R, 

~d I r 


+ 


k\ 


z ( / + 1 ) 


R, =0 


l R J 


(28) 


the solution of which is 


A(r)=V"'' + V 


IK, r 


K h 


{ e 1(1 + 1)] 


(29) 


Equating this result with that of Eq.(18), and letting B t = 0 as before, gives the 
approximation 


+'V 

ji(k\ r )~ Aj — — (30) 

a R 

Hence, using this approximation in Eq.(27), one has for the surface impedance 
at r = a R = a- S, 


z ' = -^ 
b \ 


1 


/2 


1- 


/(/+i) 

k\a\ , 


. 1/2 


(31) 


The modal value of the combination /(/+ 1) will be detennined by the behavior of the 
radial equation for region 2. 
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c. Solution for Region 2 


Consider now the scalar potential within the regolith, region 2. One thus has from 
Eq.(16) and the coordinate geometry depicted in Fig. 1 


fMh U /('+■) 

/\<2 


dr 2 


■ + 


r~ J 


R t (r) = 0, k\ = s 2 kl, a-8<r < 


(32) 


where 


' - £> + ( a “ <>) - + 1 


(33) 


Using Eq.(33) in Eq.(32) and changing variables from r to gives 


d% 2 (>) 


+ 


/(/+l) 

(^2 + ««) 2 




(34) 


Since the condition <^ 2 « a R prevails, one can follow the example given in Appendix B 
and expand the denominator of the second term within the parenthesis, viz., 




-2 


f A 
1 + ±1_ 
V 


1 - 2 — 

V a R ) 


(35) 


thus allowing Eq.(34) to be written as 


'(<+!), 2/(/+l) C 

/t 2 2 -r 3 A 


r/c : 


■ + 

2 V 


3 5 2 


A (0=0 


(36) 


Changing variables once again using 


^ 2/(/ + l) V/3 

V a* 2 


#2 


(37) 


Eq.(36) assumes the simple canonical form 

d 2 R, 

dxl 


-tR, =0 

h 


(38) 


where 


= {t 2 (l)-x 2 ) 


(39) 
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and 


*2(0= 


4 V 7,2 /(/+01 

2/ (/ + l)J r 2 al 


(40) 


The solutions to the differential equation of Eq.(38) are in terms of the Airy function 
which will be discussed much more completely in the next section. 

The parameter t 2 (I ) detennines the modes of the radial function given by Eq.(38). 
To examine the relationship between this quantity and that of the product /(/+ 1), one can 
write Eq.(39) in the following fonn 


(L-A) 3 =2 2 4(I)L 2 (41) 

where 

L = 1(1+1), A = k 2 2 a R (42) 


Solving the cubic equation that results from the expansion of Eq.(41) and selecting the 
real solution out of the three roots that correspond to the total solution gives a 
relationship that involves many terms. However, in the event that one has t 2 (l) « A 13 , 
the solution can be expanded to give 


L = A 


1 + 2 


1/3 


A 


1/3 


2 + 


144 

54A 


1/3 


t\ + - 


+ 


8 t\ 

3(2 2/3 )A 2/3 




l + 2 2/3 -^+r + — 2 1/3 q 


A 1/3 3 A 


2/3 


+ • 


(k 2 a R y 


1 + 


l k 2 a R ) 2 


, t 2 (l)« (k 2 a R f 3 


(43) 


This is a most important result that will be employed many times throughout this work; 
the combination /(/ + l) can be approximately replaced by (ka j in which the k or a 

values are peculiar to the region being considered. In some other cases, as will be 
indicated, the entire expression of Eq.(43) will be used. In fact, one can return to Eq.(3 1) 
and use this result to obtain for the surface impedance at r = a- S= a R 


f 


%=- 


4 


n 


1 _a 

. kl 


1/2 


(44) 
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Although Eq.(40) is the formal solution for the radial component of the potential 
within the regolith, one can obtain a simpler, approximate solution which exploits the 
thin layer that defines the regolith. Returning to Eq.(36) and using the fact that a R ~ a 
(i.e., S« a) and q 2 « a, one can write 


d 2 R h (r) 
d£ 


+ 


t 



1(1 + 1) 


K W=o 


(45) 


the solution of which is simply 


K (r)= A 2 e' 


I (r-a R ) 


+ B 2 e 


-iK h (r-a R ) 


a R <r < a 


(46) 


where 


f 



/(/+!) 


1/2 


(47) 


As before, the determination of the values for /(/ + 1) is connected to the radial equation 
which prevails in region 3, the exosphere. Substituting Eq.(46) into Eq.(13) yields for the 
scalar potential in this region 


U H (r,e)= ^ P, (- cos 6) 


(48) 


Once again, from Eqs.(10) and (11), the corresponding fields are 

M) 


4 (+< >**- ^ + j, Yl ( "" > jg (-cos e) 


do 


(49) 


( r ’d)= +ik o 


2 A 2 e iKl ^ aR KB 2 e iKhir ~ aR) 


dP, (- cos 9) 


do 


(50) 


These components must be matched to the previously calculated surface impedance at 
r = a R . To this end, the surface impedance within region 2 at the surface r = a R is 
defined by 


Z 2 = 




H„+.e) 


( 51 ) 


Hence, using Eqs.(49) and (50), one obtains, after simplification, 
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( 52 ) 


K l 2 (A 2 ^l) 

k(fi (A 2 + ^2} 


Since Z X =Z 2 , equating Eqs.(44) and (52) gives 


^2 B 2 _ k 2 ^ 
A 2 + B 2 k , 


12 


where 


A j2 — 


{ \l/2 / 0 Xl/2 

• k 2 

2 | | 1 _ _*2_ 


S 2 

v^i ; v 




From Eq.(53), one obtains the amplitude ratio 


B-, K, - LA 


2 12 


^2 *"/, + k 2 A l2 


(53) 


(54) 


(55) 


d. Solution for Region 3 in the Special Case of No Exospheric Plasma 

At the outset of deriving the full solution for the charged exosphere, the simplified 
case where s 2 (j ; 3 ) = 1 in region 3 will be considered first. That is, only the effects of 

diffraction and the presence of the regolith will be considered here. This will allow 
establishing of several mathematical techniques, especially those having to do with 
interfacing the source conditions into the solution, which will prove useful in the more 
complicated case in which s 3 (£,) is given by Eq.(7). Equation (17) becomes, in this 
case, 


d % 3 (r) 

dr 2 


+ 



/(/+0j 


K (/)= 0, k; = k \ , 


a < r, 4s=r — a 


(56) 


As before, one can adopt the approximation 


r 2 ={a + &)' " 


-2 


= a 




(57) 


allowing Eq.(56) to become 
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( 58 ) 


d% (r) 
dr 2 


+ 



'fi±^ + 2^)^(r)-0 

a~ a ) 


This is the same fonn as Eq.(36); thus, following the same procedure as earlier, one can 
write 

d 2 R, 

<59) 

where 

r = (^(0“ x 3) (60) 

with 


, 2/3 




21(1 + 1 ). 


^3 


1(1 + 1 )) 


(61) 


Here, 


x, = 


r 2/(/ + iy t 

3 $3 

\ a 


(62) 


A solution to the differential equation of Eq.(59) is in terms of the Airy function w 1 (t) or 
Ai(x) [11,12,13] 


i 7i ^ . 2 n ^ 

R l} (* 3 )= A 2 w l (t)= A 2 2\[ire ~*Aiyte 3 J (63) 

where w 1 (t ) is the Airy function as defined by Fock [12] and Ai (• • •) is the Airy function 
as originally defined by Miller [14]; the latter function is used in most mathematical 
references [13]. (There is a second solution to Eq.(59) denoted by w 2 (t) or Bi (• • •) but is 
ignored here due to violation of boundary conditions at t — > 00 .) 

As before, following the procedure leading to Eq.(43) with t 3 (/)« (k,a ) 3 
prevailing, one has 


/(/ + !)* (k 3 a) 2 


r 2^ 

+ 1 7 I h 

\ k 3 a) 


2/3 



(64) 
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Hence, using Eq.(47), Eq.(55) becomes 


B 2 _ (k 2 - k;J~ - k 2 A 

^2 (k 2 -kiy 2 +k 2 A 


(65) 


Now, using the general solution given by Eq.(13), one has for the scalar electric potential 
in region 3 


u. ( r,$)=A , *^ P, (-cose) 


( 66 ) 


where d 3 is another constant. But, by Eqs.(62) and (64), 


( 2kV\ 


1/3 


x 3 = 


iv) (r - a) 


P(r -a), P=k 2 y 


( 2^ 


1/3 


k 2 a 


Hence, Eq.(66) becomes 


U <r,e)=A, ^2^k^) 


P,(-COS0) 


which, when used in Eqs.(10) and (11) applied to region 3 gives, respectively, 

E e Xr,0)---p(A)„;(,,-p(r-a)Pt^ 


(67) 


( 68 ) 


(69) 


and 


H \ (a 0) = +ik (t n: h , (f 3 - p(r - 


(70) 


Rather than defining two more impedances across the surface r = a and equating, one 
can match the electric and magnetic fields across this boundary straightaway. Thus, 
equating Eqs.(49) and (69) at r = a gives 


UK, ] 

SI’ 

1 

1 bo 

1 to 

1 

Oq 

1 

V ) 

1 

hK 

to 

1 


( 71 ) 


Similary, equating Eqs.(50) and (70) at this boundary yields 
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(72) 


njA 3 w t (t 3 )= A 2 n\ 



+ ^e‘^ S 
A 2 


Finally, taking the ratio of these two results gives the transcendental equation for t 3 


w' (t 3 )~ qw,(t : ^)= 0 

where 


(73) 


q = 


J J_Yj_Yk, 


1/3 


K, 


iK, S B 2 -iK, § 

e 2 2 

A 

e iK h s + ^L e iK h 5 

A , 


* ( ,=@-*0' ! =(*2-*.T 


(74) 


Given all the parameters that define the propagation problem, one determines, with the 
help of Eq.(55), the prevailing value for q . The solution to Eq.(73) gives a spectrum of 
values for the quantity t 3 which determines the propagation modes of the 
electromagnetic radiation that can exist in region 3. It remains now to introduce a source 
of the radiation and couple it to the expression for the scalar potential, Eq.(68). 

To this end, return to Eq.(63) and nonnalize the solution at the boundary r = a . 
Thus, defining 


K fcL = ( ? 3 )= 1 


(75) 


determines the constant A } . Using this result in Eq.(68) and remembering that there are 
many values t 3 , j = 1,2,- ••, for t 3 that issue from Eq.(73), one can rewrite Eq.(68) by 
employing superposition and obtain 


> l l (^3 j 

where A. is an array of constants, one for each mode, to be determined by the source 

boundary conditions. These constants are defined to be the modal amplitudes, i.e., the 
amplitude of each of the modes supported along the lunar surface. 


U„ M)= T A J 
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i. Derivation of the Equation for the Modal Amplitudes 


In an effort to isolate the A . factors, one can attempt to treat the quotients 

Wi(f 3 . -x^jw x ^ as a set of orthogonal functions; this could be possible since the 

Airy functions come from a second-order differential equation. Therefore, consider the 
following development 


r w i 

J rU ej (r,0)— 




1 

P,(-cos6) 



(77) 


remembering that x 3 = fig . , g. = r — a . The integral occurring on the right side of 
Eq.(77) will now be reduced using all the known properties of the functions w l (•••). 
These properties are, from Eqs.(59), (60), and (73) 



(78) 


and 



dx 3 




(79) 


The first relation is the definition of w, (• • •) and the second is the associated boundary 
condition. Multiplying Eq.(78) by the function w 1 (t r -x 3 ^ and, conversely, multiplying 
Eq.(78) applied at t. by w 1 -x 3 ^ and subtracting the latter from the fonner gives 



Letting x 3 = [iz y and integrating Eq.(80) with respect to z yields, after integrating by 
parts, 



-HO 
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( 81 ) 





Hence, by Eq.(79), one has 

oo 

| Wj (t , - fiz 3 ^ (t h - fk, = 0, t 3j * t 3f (82) 

o 

Thus, the Airy functions as defined by this problem are orthogonal. For the case in which 
=t, , one must return to Eq.(8 1) and consider 

*7 °j' 


00 



Employing L’Hospital’s rule and differentiating both the numerator and denominator by 
t 3 gives 




dw\ ( 

1 3 -x,'] 
kJ k 

)dw \ ( 


dt , 

J j 


1 

ro 

•8 


1 

1 


1 * 3 , 


(t 



(84) 


Finally, using Eqs.(78) and (79), Eq.(84) reduces to 

oo 

| w i ( 3 , - As )= - R (f 3 , - r ) i 'i (f 3 , ) (85) 

0 P 

upon noting the fact that the tenn at infinity vanishes due to the asymptotic properties of 
the Wj (• • •) functions. Hence, returning to the original expression on the left side of 
Eq.(77), one has from Eqs.(82) and (85), 
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( 86 ) 


oo 00 


7=1 0 



B j, X 3 

M 

h X 0 

Wl 

i) 

Wl 

(O ' 


d^=A r 




thus isolating the modal amplitude A y , viz, using Eq.(86) in Eq.(77) gives 



(87) 


ii. Evaluation Of the Modal Amplitudes for a Given Source Distribution 

The evaluation of the expression of Eq.(87) must be done by specifying the source 
of the electromagnetic radiation in terms of its electric scalar potential. Using a 
technique that dates back to Sommerfeld [15] which is briefly detailed by Wait [16], let 
the source of radiation be a vertical electric dipole of length L be placed at 0 = 0 and 
subtending an interval g, < g : , < g^ + L along the radial coordinate. The integral in 

Eq.(87) must now be evaluated in the immediate neighborhood of this source. The 
electric scalar potential for an electrically small dipole source is, reckoned to the radial 
coordinate [17], 


tj ik 3 D 

U (r, 6) = i , D = (r 2 + r 2 - 2 r r cos &) 

630 v ’ co D v J 


( 88 ) 


where r s = a + , J is the current density within the dipole, co is the angular frequency 

of the radiation ( co = k x c ), and D is the distance to the dipole from the location r . The 
angle 0 is the angular displacement between the two radii /; and r . In the evaluation of 
the integral of Eq.(87), the point r — > r s thus requiring the limit 0 — > 0 . With these 
provisos, substituting Eq.(88) into Eq. (87) gives 



\lg, = limz 

e->o 


JL ) ]( A H; i (^3, -P(r-A)e k ' n ; 


(89) 


Since r s is taken to be the distance to the center of the dipole of total length L , the length 
the dipole subtends along the radial coordinate can be written as 

, v -LI 2 LI 2 

r = r (l + ri), <r/< (90) 

r s r s 

Equation (89) can then be written as 
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where the second result issues from the fact that ( L / 2 )/r s « 1 and cos# ~ # for # — » 0 . 
Since the dipole is taken to be electrically small, i.e., L « A, the exponential within the 
integral of Eq.(91) can be taken to be equal to unity. In this case, the remaining integral 
can be performed to give 


J rU e 3 (r,0) 


in j 

w < kh ) 



Wi (#, - p(r s - a 

V CO j 

(>, ) 

JL \ W A 

-P(r. ~ a )\ 

CO ) 

"TO 


Tim In 

0—>O 


L/2 


^ L/2 


L/2 


X 




L/2 


#^o 


(92) 


where the last result is obtained by asymptotically evaluating the limit of the logarithm. 
Therefore, Eq.(87) becomes 


A, = 


P A 


W i( 

6, -P(r* 


(s q ~ J 

^(-cos#) 

V CO J 

w A 

:o 

\ 

► 




#^0 


(93) 


However, the condition ^ — > 0 is prevailing for this expression; hence, the Legendre 
polynomial in the denominator must be evaluated in this limit. Using the well-known 
[18] asymptotic form 


^ sin In, (1-cos#) sin Ik, 

limP ; (-cos#) = In * In# 2 

/v ’ n \ 2 ) n 

the expression for the modal amplitude finally becomes 


n(3 

i | 

(JL'\ 

w A 

1 

35s 

lx" 

-a)) 

(j, -i 2 ] 

^in/;r 

V CO J 

Wi( 


1 


(94) 


( 95 ) 
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iii. The Electric Scalar Potential and Associated Radial Electric Field 


Returning to Eq.(76) for the electric scalar potential in region 3 and substituting 
into it Eq.(95) gives 




1 


. ^ 

1 

-«)) 

h~4 2 ] 

) w A 


i 


/ W 

A (t h ) Vr) 


At this point, another approximation is in order. From Eq.(64), one has 

/ « (A: 3 a) 


/ \ 2/3 

, f 2 ^ 

1 + I 

\k,a) 


— 1 1 /2 


, \ 1/3 

*3« + |^fJ *3, 


(96) 


(97) 


Thus, the value for / is such that / » 1 . In addition, the roots of Eq.(73) will be 
complex and thus, in general, / will possess a large complex part. One can then employ 
another asymptotic expansion of the Legendre function for large orders [19], viz, 


lim P, (- cos 

/| — ^ 00 ‘ V 



(98) 


The use of this expansion in Eq.(96) results in a rather unwieldy expression which would 
be difficult to treat. At this point, the following method is adopted [20]. From the 

identity sin cp = (e up - e , one has that if Im \(p } » 1 , the first exponential will 

decay away giving sin^> = /2i . Hence, since Im {/}» 1 as per Eq.(97), Eq.(98) 

can be written 


lim 73 (- cos #) = 

1/1— >00 V ' 


nl sin 0 


1/2 



^-/[(/+l/2)(;r-<9)+7r/4] 


(99) 


Similarly, 


sin In = 

2 i 

Using Eqs.(97), (99), and (100) in Eq.(96) yields, after simplification, 


( 100 ) 
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1/2 


- i 37 r /4 ik-idO 

e e 3 


Z 


r 


y=i V fc 3 


1 3 

j w i( 

-p( r s 

-«)) 

k( 

% ~P( r 

~ a i) 

-rj 

! Wl( 


) w \ ( 

:o 

I 


ix ( 0 ) t 3 . 


( 101 ) 


where x(0) = 8(kxi/2'j 3 is the scaled angular displacement from the source, sometimes 
called the reduced horizontal distance from the source [21]. 

With the electric scalar potential in hand, one can now obtain expressions for any 
component of the electromagnetic field along the lunar surface. What is of interest here, 
however, is the electric field component existing along the radial direction; this is given 
by Eq.(A40) of Appendix A, viz, 


E (r, 0) = — fsin 6^^ 

J rsmOdQV dO 


(102) 


But, by Eqs.(B18) and (B26), one has 

i df. du) 


sin# — 2- +/(/ + iVy =0 

sin# 80\ 80 J V 7 63 


(103) 


Therefore, 


E„ (r,e)= It/ (r,e)= fef u H (r,e) 


(104) 


Hence, for the radial electric field in region 3 above the lunar surface with no exosphere, 
one has 


E r (r,0)=-2i7re- ,3rr ' 4 

oo 

z 


(V) P\ wj[— Jl 


1/2 


co ) 


ik-,aO 


{ * A 

J w i( 


- a i) 

k( 

/.-A'- 

-a)) 

U ,-q 2 ) 

Wl( 


I w A 

:o 

t 

1 


ix { 6 )t 3 


(105) 


Finally, letting the height to the source from the lunar surface be =r s - a and 
similarly the height from the surface to the receiving point be <^ 3 = r - a and noting that 

r 2 = + a J ~ a 2 , one has, upon using the definition for J3 and rearranging factors 
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E.ik-ey E Mv(x(0U,,,4,,‘i) 


(106) 


where 


V 


co 

(x{6)^ x ,c x q)= le n! * (xx(0)J~ Yj 


1 ' 


f 

t 3 

-A,.) 

k( 

f 

t 3 

-K,) 

-rj 

1 Wi( 

' 5 
t, 

) 

) H l ( 


i 


ix(0)t 3 . 


(107) 


is the ‘attenuation factor’ [22] and 


E 0 (O)=ik* 


JL\ e ihae 

a> 7 a(#sin#) 12 


(108) 


is the non-diffracted electric field strength. Equation (107) is the celebrated expression, 
originally obtained using other methods by Fock [23], for the diffracted electric field 
strength of an electromagnetic field along the surface of, in this case, the moon. The 
effect of the regolith is convolved in the parameter q given by Eq.(74) as well as in the 
solutions t 3 of the modal equation, Eq.(73). It remains now to solve Eq.(73) for the 
spectrum of t 3 values and incorporate them into Eq.(107) for evaluation. 


iv. Numerical Evaluation of the Modal Equation and of the Electric Field 

One must necessarily commence with the numerical evaluation of Eq.(73) in 
order to obtain the array of t 3 values to be used in Eq.(107). It is useful to decompose 

these values into their real and imaginary parts, viz., 

t 3j =Xj + iyj (108) 

Since the evaluation of these expressions will be done using Mathematica, and 
Mathematica only uses the Airy functions as well as its derivative At '(•••), one 

must use the representation given in Eq.(63) and rewrite Eq.(73) as 

e 3 Ai' j^(x y + iVj )e 3 j - qAi ^ (x y + iy } )e' 3 j = 0 (1 09) 

By plotting zeros of the real and imaginary parts of this expression with respect to x . and 
y , , one can immediately identify its roots to be at the points where the zero contours 
intersect one another. 

For example, take the frequency of operation to be 100 KHz; the relevant 
electrical properties are taken as displayed in Table. 1. Using these values in Eqs.(54), 
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(65) and (74) (noting that k 3 = k 0 here since e 1 = 1 as the exosphere is absent), one 
obtains that in the case of no regolith ( 5 = 0 ) 

q = -0.857 - 0.936/ (110) 

Using this value in Eq.(109) and plotting the zeros of the real and imaginary parts of this 
expression gives the plot displayed in Figure 5. 


100 KHz 

q = -0.857-0.936; 



Figure 5. Real (x) and Imaginary (y ) Parts of the Argument of Eq.(109). 

Here, eight modal values of t 3 are displayed. The first root A is the dominant mode 
since it has the smallest imaginary part and will thus yield the smallest attenuation 
occurring in the e' A ^ 3 ' factor of Eq.( 106). Due to this attenuation contributed by 
hn -jf 3 |> 0 , only the first four or five values of t 3 need be used in the evaluation of 

Eq.(107) as the factor will quickly decay away for large hn <jt 3 j All other 

calculations similarly follow. In some cases, one will encounter the first root in which 
hn j/ 3| j< 0 is an unphysical result and is thus an extraneous root. Hence, the root must 

be rejected and the next root above it for which hn |> 0 becomes the root 

corresponding to the dominant mode. 

In the results that are to follow, transmitter and receiver antennas are, 
respectively, 20 m and 2 m above the lunar surface ( = 2000 cm, q, = 200 cm), and 

the transmitter is taken to have a power output of 1 W. The transmitting dipole length is 
10 m. Cases for the transmitting frequency of 50 KHz, 100 KHz, 200 KHz, 500 KHz and 
1 MHz are shown in Figures 6-10 which show electric field strength in volts/meter versus 


NASA/TM— 2008-2 1 5463 


26 


distance from the source in kilometers. At 50 KHz, there does not appear a trapped 
surface wave phenomenon; each case is not too different from the others. Figure 7 
displays the rather interesting situation at 100 KHz in which the 8 - 100 m case does 
much worst than that where 8 = 0 ; here, the electric field intensity is almost an order of 
magnitude down from that of 8 = 0 at a distance of 100 Km. The regolith of 8 = 1 00 m 
attenuates the wave along the surface as compared to the other cases. However, for 
8 = 50 m, a trapped wave phenomenon seems to occur allowing the field strength to be 
more than an order of magnitude greater than for the 8 = 0 case. 


V 

mtr 


50 KHz 



<5 = 100m 


-<5 = 50m 


<5=0 


Figure 6. Electric Field Strength vs. Distance Along Lunar Surface at 50 KHz 


V 

mtr 


100 KHz 



<5= 100m 


Figure 7. Electric Field Strength vs. Distance Along Lunar Surface at 100 KHz 
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V 


200 KHz 


mtr 



Figure 8. Electric Field Strength vs. Distance Along Lunar Surface at 200 KHz 


__ 500 KHz 

mtr 



Figure 9. Electric Field Strength vs. Distance Along Lunar Surface at 500 KHz 



Figure 10. Electric Field Strength vs. Distance Along Lunar Surface at 1 MHz 


At 200 KHz as shown in Fig. 8, the trapped surface wave ceases to exist and both 
regolith cases do worse than that for S = 0 . From this example on, the regolith attenuates 
the propagating wave field along the surface as shown in Figs. 9 and 10. Hence, the use 
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of frequencies below 100 KHz would not be warranted in order to exploit the trapped 
surface, ‘whispering gallery’ wave phenomenon. However, what is most disconcerting is 
the fact that in the real case of the lunar regolith, one cannot expect it to remain at a given 
depth. In fact, the regolith depth will be some highly random function of direction and 
distance from the transmitter site. Hence, one can expect that, given the enhancement 
and attenuation of the varying regolith thickness, the electric field strength will most 
likely be approximately at the value realized with the case 8 = 0 at 100 KHz. This, of 
course, requires a more detailed investigation. 

The behavior of the signal field strength as a function of both frequency and 
regolith depth also requires a more careful investigation than given here. The analysis of 
this spherical dielectric waveguide problem would not be too different than that used to 
study the cavity fonned by the Earth and ionosphere [24]. 


e. Solution for Region 3 for the Case of an Exospheric Plasma 

The solution of the general problem defined by Eq.(17) using the permittivity 
profde given by Eq.(7) will now be addressed using the mathematical methodology 
established in Section Ill.d. Thus, one must now deal with the differential equation 


d fe) , { f „ , „ f & ~&o ) ] lQ_+j) 

( a + &J 


dg 




A(&)=° 


(no 


obtained after changing variables to the local coordinate g, . Using the approximation of 
Eq.(57) as well as the change of variables noted by Eq.(62), Eq.( 111) becomes 


d R t (x 3 ) 


dx x 


+pQ) 


k r - 

^0 c 30 


/(/ + !)' 


K (u) 


Pif)K E 3 


+ 


( x 3 — x 30 ^ 


v 30 


+ X? 


R<X x i)= 0 ( 112 > 


where p(j)= (a 3 /(2l(l+ l))J . Defining 


k 2 = k 2 r k 2 = k 2 p 

^30 — ^0^30? /V 3A “ ,t 0 d> 3A 




i-- _ 
^30 


/(/ + !)" 


fXx 3 j)= p{l)k 2 ^y 


( x x - x 30 ^ 


^30 


+ X. 


(113) 
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Eq.(l 12) can be written as 


^^+(/ s (v)- >Mh 69= 0 ( 11 4 ) 

Again, as was done in the procedure leading to Eq.(43), one needs to find the solution 
/(/ + l) to the equation 


,,, J a 3 t 

{llfl + l} 


2/3 


7(/ + l) 


-k: 


2 30 


Hence, one has 


/(/ + 1)~ (^ 30 a) 2 


/ \ 2/3 

1 + 1 T J ^ (0 


^3 ( 0 « (^" 30 a ) 2 ^ 


Using the further approximation /(/+!)« (k Vl a J , one can write 



The function / 3 (x 3 ,/) now becomes 


L (x 3 ,l)= f 3 (x 3 )= 

V Un 


+ x 


3 ’ 


K = 



Equation (114) thus becomes 


d 2 R, 3 (u) | 

dx 3 




(115) 


(116) 


(117) 


(118) 


(119) 


Although Eq.(l 19) can be treated by ‘completing the square’ of the x 3 term, it 
proves to be more expeditious to perform a Taylor series expansion [25] about the 
minimum of /) (x 3 ) . To this end, one has 



( 120 ) 
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The minimum x 3m of f } (x 3 ) occurs at f' (x 3m ) = 0 ; hence, the first relation of Eq.(120) 
gives 


X 3m X 30 I ^ 


-^30 

2 K 


Therefore, 


f3( X 3 m )= X 30{l-^ 


Thus, one has for the function f 3 (x 3 ) 


M X 3)=M X 3m) + \f" ( X 3mX X 3- X 3mJ 


X 30 I 1 


~^30 

4 K 


+ \i X 3- X 3m)~ f 3 


which gives, upon substitution into Eq.(l 14) and changing variables, 


d - R ^ + - my E)= o 

where z 3 =x 3 - x 3m and 

Now, following Fock [25], consider the change of variables 

C=fy"Je-“ H z, 

Equation (124) finally takes the fonn 



where 


(121) 


( 122 ) 


(123) 


(124) 


(125) 


(126) 


( 127 ) 
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v+y 2 =~i 


(128) 


The solutions of the differential equation of Eq.(127) are parabolic cylinder 
functions D v of order v [26, 27], viz., 


R,X()=A 1 D_ l _ v (iC)+B 1 DXC) 


(129) 


The behavior of this solution as £ — > oo renders the constant A 3 = 0 . From Eqs.(126) and 
(120) and the definition of z 3 , 


C=a(fi(r-a)-x 3m ), a = ^ 


i AK T — * * 

I — - <? , p = k 3 


rj_v 

'l k 3n a) 


(130) 


Thus, one has, from Eqs.(13) and (129), the scalar electric potential within the exosphere 


TT / Q \ A „ D v(a(P(r-a)- x 3m )) 

(r, e) = AB 3 ^ *P t (- cos 0) 


(131) 


with v given by the solution of Eq.(128). The associated electric and magnetic fields are 
given by, using Eqs.(10) and (1 1), 

E e 3 M)= a PD[. (a (p(r ~a)-x 3m ) f P ‘ ^° S ^ (132) 


H +, ( r ^)= iK n lAB 3 {^\ D v (a(p(r- a)- x 3m ) f P ‘ ^ (133) 


Following the same process as before and equating these components at the boundary 
r = a with those in region 2, using Eqs.(49) and (50), gives 


AB 3 a[iD\{-ax 3m )= i/c,A 2 e' Kh - S - MV? 


(134) 


An] B 3 D v (-ax 3m ) = A 2 n] e' K ' 2 ° + y- e' K ' 2 ° 


(135) 
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Here, remembering Eqs.(47), (55), and (116), 


and 



( 136 ) 


A ( 

7 2 7.2 

n 2 _ ^30 y 

^ /2 

) _ ^- 2^12 

7, ( 

k 2 - k 2 

K 2 a; 30 

a /2 

) ^ 2^12 


(137) 


with A 12 given by Eq.(54). Taking the ratio of Eqs. (134) and (135) gives the modal 
equation 


D l (~ax 3m )- qD v (~ax 3m ) = 0 


(138) 


where 



The solution of Eq.(138) for the propagation modes within the exosphere is not 
with respect to the arguments -ocx 3m for they are specified by the problem; the solution 
is with respect to the order v . These spectra of values i / . are then used in Eqs. ( 128) and 
(125) to determine the corresponding values of t 3 for use with Eq.(l 16). 

Normalizing the radial solution given by Eq.(129) (with A 3 = 0 ) at the boundary 
r = a yields 


R h (a) = B 3 D v (~ax 3m )= 1 (140) 

Hence, the full solution for the electric scalar potential is, using Eq.(13 1), and noting that 
one is dealing with an array of modal values v. 


v„{r,e)= 




7 = 1 


D v. (~ ax 3 m) 


- E’ / (-cos^) 


(141) 


As before, the values of the amplitudes Aj of the modes supported along the lunar 

surface will be determined by the boundary value problem defined by the source of the 
electromagnetic radiation within the exosphere. 
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i. Derivation of the Equation for the Modal Amplitudes Within the 
Exosphere 

The procedure to obtain the values of the modal amplitudes will follow along the 
same lines as established in Section IILd.i above. Thus, one attempts to decouple the 
amplitudes A j from the sum of Eq.(141) by forming the integral 


\rU e Xr,0) 

0 


D Vr ( a (ff£»-*3m)) 1 

D v f (r aX 3 m) 3 P,(-COS0) 

_y A r D v, - x im ))D Vj (a (pc, - x 3 , 

j = 1 0 D v r (r aX 3 m) D v (- a X 3; „) 



(142) 


and attempt to ascertain the orthogonality of the parabolic cylinder functions appearing 
on the right side using the established properties of these functions dictated by the 
problem at hand; that is, the defining equation of Eq.(127), 


■rty (£) 

d£ 2 + 



(143) 


and the condition of Eq.(138) 


= 0 (144) 

g=-ax 3m 

Thus, the procedure established in IILd.i will be followed; multiplying Eq.(143) by 
D v (C) and subtracting from this result Eq.(143) for v y multiplied by D v (^), changing 

variables using C, = a/?z 3 - ax 3m , integrating with respect to z 3 , and performing an 
integration by parts on the portion with the derivatives yields 

oo 

= ( v r~ v j 

z 3 =0 

Using Eq.(144), this gives 


)jD rr (Cp,.(0dz, (145) 




dD (£) dD Vi (Cf\ 

VjKJ -d Vj (c) 


dC 


dC 


dj\ fe) 

dC 


~9D y (d 


\D Vf (Cp Vj (C)d Z 3 = 0, Vj*Vj. (146) 

o 
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In the event that v. = v jt , one has from Eq.(145) upon applying L’Hospital’s rule 


} , . , . 1 ( dD (C)dD v (£) d 2 D 

lim D (Cp y ( C)dz 3 = — lim — J —~ D v (£) 


( 147 ) 


z-, =0 


The connection between £ and v y is given through y = v ]f + 1 / 2 - (1 / 4)p ; hence, 
%r-0<y d A r . Using this in Eq.(147) and employing Eqs.(143) and (144) 
finally gives 


I A <&> = - («X ( 7 + u h) b + X - (Kk H <->) 


f 1 V 2 ^ 




- ( V J + K)“ ^ )°v, (-«U m ) (148) 


Hence, using the results of Eqs.(146) and (148) in Eq.(142) gives for the modal 
amplitudes 


A 


j 



1 

P t (- cos 0) 


J rU e 3 ( r,0 ) 
0 


D,, (aj/A-x,, 

D Vj (~<*X 3m ) 



(149) 


ii. Evaluation of the Modal Amplitudes for a Given Source Distribution 

The source of the electromagnetic radiation in the lunar exosphere is again taken 
to be a vertical electric dipole, the electric scalar potential of which is given by Eq.(88). 
Applying this relation to that of Eq.(149) and following the procedure given in Section 
IILd.ii, one has for 6* — » 0 , 


\rU e 3 M) 


D v (a(fA-x, m )) 


D Vj (~ ax 3 m) 


d4 . 3 


jy] D A(p('-.- a )- x s ,% |n ^ s 

a) D Vj (~ax 3m ) ^ ' 

(150) 


Substituting Eq.(150) into Eq.(149) and using Eq.(94) gives for the modal amplitudes 
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i{ na 2 0x 3m ) (J_£\ ( I 1 

2 l sin nl Aco){(y 4 JaxJ-(v J + / 2 y q 2 ) D v (-ax 3m ) 


( 151 ) 


iii. The Electric Scalar Potential and Associated Radial Electric Field 

Substituting Eq.(15 1) into Eq.(141) gives for the scalar electric potential in the 
lunar exosphere 


u*(r,e) 


-_L{ na P X lm 1 f I 1 . 

2 l sin nl J l CO ) j=x [ Qfyax iM J - (v. + J/Q- g 2 J 

D Vj (<*(p(r s - a)- x 3m ))// (a(f:i(r - a)- x 3m )Y A 
D vj (~ax 3m ) D Vi (-ax 3m ) 


(152) 


Using Eqs.(97), (99), and (100) in Eq.(152), and employing the definitions for a and /? 
yields 


TT A 0\ A — />/4i f ^ | jf^l/2 ^ f f ^ ^ 1 * 30 oS 


y=i 


r 

i 

i 


5 

^ { 

[a(f:)(r s - a)- x 3m )^ 


(a(/?(r-a)-x 3m )) 

1(h) 

(«0 2 -{ 

r , 1/3 

I-? 2 ; 


D vj 






(153) 


where, from Eqs.(125) and (128), 


t 3 = / 


,.(4£T 


v 30 


'(vX)*- 


^30 


(154) 


The associated electric field for this scalar potential is given by Eq.(104); 
substituting Eq.(153) into Eq.(104) yields for the radial electric field within the exosphere 


r 2 Y 3 x^C JLV 1 ^ 


{r^Te-X^y - j 


, k 30 a sin 6 


1P ( l x2 

e ik io“0\ _ 

v r 


r 

i 

i 


5 

A, { 

[a(fi(r s ~ a)- x 3m y 

K< 

(a(/?(r-a)-x 3m )) 

1(h) 

(«aJ - 1 

6 1/3 

i/ + y 

K> 7 2) 

I-? 2 ; 


D v, (~ax 3m ) 


D Vj(r aX 3») 


*(^3 


( 155 ) 
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Finally, letting g 3 be the distance from the lunar surface to the point of reception, i.e., 

= r- a , and r~ 2 = +aj ~ a 2 with =r s - a , Eq. (155) gives, after 
rearrangement of factors and simplification 

E ri (g,e)=E„(g)v(x(g),g,gq) (i56) 

where the attenuation factor is 


v(x(9), 4 , £ q y -2ie M4 K ]n (xx(e)f ( ^ 

V x 3n 


j = i 


r 

i 

i 



1 ^ I 

f (f%3, 

K < 

«(/^ 3 , - )) 

1(h) 

( ax ,J - 1 

r i/> 

y + y 

K> /2) 

)-q J 


( ccx 3m ) 


(~ aX 3m) 


(157) 


and the non-diffractive field strength is 


E 0 (9)=ik 2 30 


JL \ e ik ™ ae 

0) 7 a (9 sin 6 1 ) 1 " 


(158) 


Equations (156) -(158) describe radio wave propagation within the lunar exosphere and 
takes into account diffraction around the surface of the moon as well as the effect of the 
regolith. 


iv. Numerical Evaluation of the Modal Equation and of the Electric Field 

At the outset, one must evaluate the spectrum of roots v j that satisfy Eq.(138). 

To this end, it is required to provide a method for the evaluation of the derivative 

D' (• • •) since Mathematica does not supply a routine as it does for the derivative of the 
Airy functions. To circumvent this, one can use recursion relationships that prevail for 
the parabolic cylinder functions [28], viz., 

D[{x) + (y 2 yD v {x)- vZV, (*)= o, D' v {x)-(y i )cD v (x)+D v+l (x)= 0 (159) 

Adding these two relationships, solving for D[,(x ) and using this result in Eq.(138) gives 
for the modal equation, 

vD v _, (-ax 3m )- D v+1 (~ax 3m )- 2 qD v (~ax 3m ) = 0 (160) 
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which must be solved for v . As done before, the contours of the real and imaginary parts 
of this equation will be plotted with respect to the parameters cr and z where 

v j = a j +iv j (161) 

are the roots of Eq.(160) where the two contours intersect. 

As an example, consider an operating frequency of 500 KHz (sufficiently above 
the 220 KHz cutoff frequency as shown in Fig. 3.). Taking a the regolith thickness of 
8 = 50 m, and using the parameters discussed in Section II. b, one finds, upon 
implementing Eqs.(8), (54), (137) & (139) 

q = -3.289 - 0.210/ (162) 

Using this and the corresponding values of a and x 3m in Eq.(160) yields the following 
contour plot for the roots v . of the equation as shown in Fig. 1 1 . The first root is noted 
to occur at v l = -0.637 - 0.658/ . Here, one must be careful that the roots v. chosen are 

r i ix(o)t 3 

such that, using Eq.(154), Im <f 3 >> 0 so as to have a solution that decays with e ,J . 

In this case, v 1 is such that Im |< 0 . Hence, this root is an extraneous one and must 
be neglected. The second root v 2 = -0.357 + 0.378/ satisfies the corresponding 
hn{/ 32 } > 0 and thus becomes the first root of the mode spectrum i /. and is the dominant 
mode. 


500 KHz 

q = —3.289 — 0.210/ 



Figure 11. Real (cr) and Imaginary (r) Parts of the Orders v of Eq.(138). 


As before, assume a transmitter height above the lunar surface to be 20 m and that 
of the receiver to be 2 m; the transmitter output power into the 10 m length dipole is 1 W. 
Taking the three regolith thicknesses of 8 = 0 m, 8 = 50 m, and 8 = 100 m, and using 
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Eqs.(139), (160), (154), and (156) yields the results shown in Fig. 12. Here again, the 
electric field strength in volts/meter is plotted versus propagation distance in kilometers 
along the lunar surface. The deleterious effect of the regolith, similar to the case of Fig. 
9, is seen. 


500 KHz in Exosphere 

mtr 



Figure 12. Electric Field Strength vs. Distance Along Lunar Surface at 500 KHz 


The behavior of the effect of the regolith is the same as it is for the case of no exospheric 
plasma save for the fact that the field strengths are more than an order of magnitude 
greater for the exospheric case. The 1 MHz example shown in Fig 13 again follows that 
of Fig. 10 for the case of no exospheric plasma. Thus, the existence of the regolith is 
deleterious to radio wave propagation along the lunar surface in the presence of an 
exosphere. 


1 MHz in Exosphere 

mtr 



Figure 13. Electric Field Strength vs. Distance Along Lunar Surface at 1 MHz 
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V 


2 MHz in Exosphere 


mtr 

0.001 
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10- 6 
10 - 7 

0 20 40 60 80 100 

Figure 14. Electric Field Strength vs. Distance Along Lunar Surface at 2 MHz 



— 3 MHz in Exosphere 

mtr 



Figure 15. Electric Field Strength vs. Distance Along Lunar Surface at 3 MHz 


For 2 MHz shown in Fig. 14, all regolith examples perform about the same 
whereas for 3 MHz shown in Fig. 15, the cases again separate in which the regolith 
attenuates the wave as compared to the case of 8 = 0 . 


IV. Conclusions and Discussion 

It has been shown in the examples considered in this study that the presence of a 
lunar regolith is deleterious to transhorizon radio wave propagation above frequencies of 
100 KHz. However, when an exospheric plasma is present, the cutoff frequency for 
communications is ~ 200 KHz, below which electromagnetic waves do not propagate 
through the plasma. Thus, unless the exospheric plasma is nonexistent, the lunar regolith 
provides no advantages for surface wave propagation. However, the exospheric plasma 
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as assumed here does actively support long-range transhorizon communications for 
frequencies from 500 KHz to 2 MHz, the lower frequencies being favored over the larger 
ones. Given the fact that the 500 KHz frequency is close to that of the exospheric plasma 
frequency, and allowing for variations of the plasma properties as assumed in Section 
II. d, it would be a reliable choice to employ communications frequencies ~1 MHz for 
lunar transhorizon communications. 

These conclusions, however, are based on an exosphere model composed 
essentially of one data point, i.e., the data of Fig. 2. It is clear that much more 
information of the exosphere is needed before reliable communications system 
predictions can be made. For example, the variation of electron concentration with 
respect to height linked to the spatial variation across the dayside surface of the moon is 
required. Some very early considerations of these points can be found in the literature 
[29]. The exosphere model of Section Il.b and the propagation model of Section like are 
general enough to accommodate a wide range of exosphere scenarios. Additionally, 
information of the lunar regolith in terms of its spatial variation with respect to 
displacements along the lunar surface as well as its electrical properties and that of the 
core material is needed to assess the impact on surface wave propagation. Finally, it 
must be remembered that the examples given here are for the worst case of an electrically 
small vertical electric dipole operating at 1 W. The use of a directional transmitting 
antenna and larger power will, of course, substantially improve performance. 

Although there is much data to be obtained, this study demonstrates the potential 
for low frequency, long-range transhorizon communication along the lunar surface and 
provides an analysis procedure for its assessment. It hopefully will provide impetus for 
further investigation and consideration for low data rate communications for contingency 
purposes during lunar exploration activities. 
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Appendix A 


Derivation of the Electromagnetic Field Components in Terms of A Scalar Potential 

for TM (Electric) Waves 

By way of a brief review, the components of the electric and magnetic field 
vectors of a propagating wave field along a spherical planet or moon will each be related 
to a single scalar potential. The particular case for transverse-magnetic (TM) waves, or 
electric waves, will be considered. The equations derived in this Appendix will fonn the 
basis of analysis of electromagnetic wave propagation in a spherically layered 
environment. 

The Maxwell Equations in traditional cgs units are given by 


- - 1 ED Ax - 

V xH = — / 

c dt c 

(Al) 

V x E + — — = 0 
c dt 

(A2) 

V • D = Axp 

(A3) 

VB = 0 

(A4) 


which relate the electric and magnetic field intensities, E and B, to the electric and 
magnetic inductions, D and B, in the presence of an external current j and charge 
density p. Associated with these are the constitutive relationships involving the 
prevailing electric permittivity a, the magnetic permeability p, and the conductivity a, 
viz, 


D = sE, B = pH, j = oE (A5) 

In the propagation situations considered here, one has p= 0 and p = 1. Using Eqs.(A5) 
in Eqs.(Al)-(A4) with these values gives 


- - scE Ax - 

V x H = — oE 

(A6) 

c dt c 

VxE + -— = 0 

(A7) 

c dt 

V ■ E = 0 

(A8) 

VH = 0 

(A9) 
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Furthermore, the field components E(r,t ) and H(r,t) are taken to be time hannonic, 


E (r, t) = E' (r )exp (-zmt) 

(A10) 

B(f,t)= 2?'(F)exp(-z<z>t) 

(All) 

Substituting Eqs.(AlO) and (A1 1) in Eqs.(A6)-(A9) yields 


- co - An - 

VxH' + i — sE' = — <jE' 

(A12) 

c c 


V xE' -i — H' = 0 
c 

(A 13) 

V-E' = 0 

(A 14) 

V • H' = 0 

(A 15) 


In this case, Eqs.(A14) and (A15) are essentially the result from Eqs.(A12) and (A13); 
hence, the former will be omitted from further treatment. Thus, the working relationships 
for the electric and magnetic fields are 


Vxtf' = f — o-i — ^lE' (A 16) 

vc c ) 

VxE' = +i—H' (A 17) 

c 

In order to capture the essence of the parameters in these equations and their impact on 
wave propagation, it is beneficial to obtain the wave equation from Eqs.(A16) and (A 17); 
in the usual way, one obtains from Eq.(A17) 

VxVxE' = v(v-E')- V 2 E' = +z — V x H' (A18) 

v J c 

Using Eqs.(A14) and (A16) and rearranging terms results in 

V 2 E' + k 2 E' = 0 (A 19) 

where 

( 4 na\ V1 

k = k 0 \ £ + i = k 0 n (A20) 

v co ) 
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Here. k 0 = <x>/c is the free-space wave number and n is the refractive index of the 
propagation medium. 

Given these considerations, Eqs.(A16) and (A 17) can be written 


VxH' = -ik 0 n 2 E' 

(A21) 

V x E' = +ik 0 H' 

(A22) 


It is from these two relations that one can obtain individual expressions for each of the six 
components comprising the electric and magnetic field vector components solely in terms 
of two scalar functions, the electric and magnetic potentials. To this end, the solution in 
tenns of the electric potential for a problem in spherical coordinates will be found. In 
spherical coordinates, 


VxH' = 


/ j N 


r sin 6 A [36 


^(rsin 0H t y-^(rH e ^ 


r — 







(A23) 


where r, 6, and <f) are the unit vectors in the radial, colatitude, and azimuthal directions, 
respectively, and H r , H e , and H tj) are the associated magnetic field components. Hence, 
the components of Eq.(A21) are 


-ik n nE„ = 


1 


r sin# 


C (sin 0H,)--^H e 


36 


(A24) 


-ik 0 n E 0 = 


1 


rsin# 


3H„ 3 


3(j) 3r 


(rsin OH^) 


(A25) 


- ik 0 n z E ^ = I - 
\r 


3(r He y 0 JL 

3r y e) 36 


(A26) 


Similarly, for Eq.(A22), one obtains 


-ik 0 H r = 


1 


rsin 6 


3E„ 3 / . 


36 


(sin 6E^) 


{All) 


ik 0 H e 


1 


rsin 6 


3 ( . n ^\ 3E 

— (rsvndE. ) - 

3r y ^ 3(j) 


(A28) 
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-%o H + 



dE r 

DO 



(A29) 


The object now is to essentially decouple these relations to obtain explicit equations for 
each of the separate field components. The process commences with limiting the analysis 
to fields of transverse magnetic type (i.e., electric waves) in which H r = 0 and E r ^ 0. 

(A similar procedure can be done for the transverse electric case in which E r = 0 and 
H r * 0 but this is not of interest in this work.) 

With H r = 0 , Eqs.(A25)-(A27) simplify to 


-ik 0 n 2 E g = - 

(;) f to) 

(A30) 

= ( 

i)Jto) 

(A31) 

II 

-STS 

*8 ^ 

(A3 2) 


The only way Eq.(A32) is satisfied is if there exists a scalar function ® such that 


sin OE^ = — , E ff = — 

* df 6 do 

Elsing the second relation of Eq.(33) in Eq.(30) yields 



(A33) 


(A34) 


Multiplying Eq.(30) by sin 0 and using the first relation of Eq.(33) gives 



sin# d 

r dr 



) 


(A3 5) 


Now assume there exists another scalar function 


y/ such 


that ® = 


1 d\j/ 
r dr 


Eqs.(A34) and (A35) to become 


This allows 


H„ = + 


ik 0 n 2 dy/ 
r d6 


(A36) 


and 
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ik 0 n 2 3y/ 
rsin# dtj> 


(A3 7) 


H a = - 


Finally, it is convenient to rescale the function y/ by 1 Jr, i.e., consider the scalar U e = — . 

r 

Equations (A3 6) and (A3 7) then simplify to 


Hi = +ik 0 n' 


36 


(A3 8) 


and 


ik 0 n 2 3U e 
sin 6 3<j> 


(A39) 


These two magnetic field components are now given in terms of the electric potential, 

U e . Given these, they can now be used in Eqs.(A24), (A30), and (A3 1) to obtain similar 
expressions for E r , E g , and E . : 


E= — 


1 


rsin# 


3 


do 


sin# 


36) 


+ - 


3 


1 3U l 
sin 6 3(f> ) 


(A40) 


E a = 


13_' 

r 3r\ 


36 


E„ = 


r i ^ 

3 

L*'.) 

V rsin #y 

3r 

l cty J 


(A41) 

(A42) 


Finally, an expression needs to be found that yields the electric scalar potential 
U . For this, Eq.(A28) is used. First, from the first relation of Eq.(33), one can write 


sin ££= — = — 
* 3(f> 3(f> 


'1' 


\rJ 




(A43) 


Using this and Eqs.(A38)-(A40) in Eq.(A28) gives, after some simplification, 

-UUfsi 

sin 0) 36 v 




3<j) r 3<j) 3r 


~ i . a 3U e ] 3 

, , ,, — sin# + — 

V r smd) 36 v 36 


r 1 3U, A 


( sin 6 3(j) ) 


(A44) 


Integrating this expression over (f) and using the fact that k = k 0 n finally gives 


NASA/TM— 2008-2 1 5463 


47 



r dr 


1 


'd f 


V r sin #7 


SO 


sin# 


do) 


1 


v #u A 


V r sin" 6\ d<f> 


+ k 2 U= 0 


(A45) 


J 


This is the wave equation in spherical coordinates for the scalar potential of TM waves. 
The solution of this equation with the appropriate boundary conditions yields, from 
Eqs.(A38) - (A42), the components of the E and H fields in the case where H r = 0 . 

Finally, through the use of Eq.(A45), Eq.(A40) for E r can be simplified. 
Multiplying Eq.(A40) by 1 Jr, adding it to Eq.(45) and solving for E r gives 


E ,=§i(’-U,)+r(rU,) 


(A46) 


The boundary conditions for the potential U e on a spherical surface follow from 
the requirement of continuity of the tangential components of the fields. For a spherical 
layer of radius a separating two electrically distinct regions 1 and 2, one has for the 
tangential component of the electric field must be continuous across the interface, i.e., 




(A47) 


Using Eq.(A41), this requires 


^ U e t 00 ) 


3- 


_ ^ U eX r )) [ 

dr 


(A48) 


i.e., the quantity 
magnetic field 


dr 


must be continuous across the interface. Similarly, for the 


H A'to H d-L 

From Eq.(A38), 


(A49) 


(A50) 


across the interface. 
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Appendix B 


Solution of the Wave Equation for the Scalar Potential in Spherical Coordinates 

In this Appendix, the well-known solution of Eq.(A45) will be reviewed for 
completeness. Hence, solutions and boundary conditions for the potential U e will be 
found from the equation 


r or 


1 


r sin 0) 


8 


do 


sin# 


89) 


+ 


1 


V 


8 1 2 u„ 


V r sin - 9 A 8(j) 


+ k~U= 0 


(Bl) 


J 


To this end, it is necessary to isolate terms with separate radial and angular dependencies. 
Multiplying Eq.(Bl) through by r 2 gives 




r -^(r u e )+r 2 k 2 U e + 


vsin6v 


d 


89 


sin# 


89) 


+ 


1 


V 8?U ^ 


vsin" 9)\ 8(j) 


= 0 


(B2) 


The form that this relation assumes suggests that one consider a solution of the 
multiplicative fonn 


U e =U e (r,0,f) = R(ry(0,f) 


(B3) 


Substituting Eq.(B3) into Eq.(B2) and dividing through by the same product R(r)Y(9,(f)) 
yields 


—r-^-r-irR^+r 2 ^ + 
R 8r 2V J 



_J_8_ 
sin 9 89 



+ 


1 3 8 2 Y 
sin 2 0/ 8(jr 


0 


(B4) 


The first two terms on the left are only functions of the radial component r and the 
remaining terms are functions only of the angular component 9. Since there is no 
functional dependency between both these groups of tenns, they both must be equal to a 
constant which is an additive inverse for each group (since they both give a final result of 
zero). That is, letting X be a constant (not the wavelength), one has 


1 

— r — t 
R 8r 2 


(ri?)+ r 2 k 2 = X 


(B5) 


iir^( sin ,Ai + 

Y) sin 9 80\ 89) 


1 3 8 2 Y 
sin 2 9) 8(jr 


-X 


(B6) 


These relations reduce to 
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-4-0)+ 

r or 


r 


e-3 

v r J 


R = 0 


1 d 


sin# 3d 


• a 3 * 
sin d — , 

3d) 


+ 


i 


vsin 0) 




3<f>~ 


+ AY = 0 


(B7) 

(B8) 


Equation (B7) gives a differential equation governing the radial function R(r) of the 
solution. Equation (B8) can be reduced further following the same procedure as 
employed above. Multiplying Eq.(B8) by sin 2 0 gives the separable relation 


• a 3 
sin 0 — 

36 


. 3Y 
sin 6 — 

36) 


+ 


rf-y 

AY sin 2 6=0 
3<j> 


(B9) 


Assuming that one can write the solution for this equation in the form 

Y = Y(6, <j>) = Q{6}Y((/)) (BIO) 


and using this in Eq.(B9) gives 


f r 

i . „ d\ 

{ . n 3&) 


sin 6 — 

sin 6 — 

V®) 

1 36 

l 3d) 


+ Asm 6 = v 


1-1 

vOy 


rfO 

~30 2 


= -v 


(Bll) 

(B12) 


These relations reduce to the governing differential equations for the colatitide function 
0(6) of the solution and the azimuthal function O(^), viz., 


1 


sin 6 36 v 


0 )m0^) + \X 

3d) 



0 = 0 


(B13) 


^0 

— ^+nO = ° (B14) 

3(p 


The analysis of Eqs.(B7), (B13), and (B14), the solutions of which collectively 
give the composite solution of Eqs.(B3) and (BIO), commences with Eq.(B14). Its 
solution is simply given by 

®(0) = A exp B cxp(-i^fv(f)^ v y 0 (B 1 5) 

0 (0) = A + B<f>, v = 0 (B16) 
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Since the propagation problem addressed here does not involve the azimuthal coordinate 
</), v = 0 ; additionally, since no particular value for (f> is specified, B = 0, and thus 


® (0)= A = const. 


(B17) 


Equation (B13) now becomes 


1 6 ^sin£— I +/L0 = O 


sin 0 d6\ d6 


(B 18 ) 


Changing variables to w = - cos 6 , Eq.(B18) takes the form 

+ AP(w) = 0 


_d_ 

dw 


V ’ dw 


(B19) 


where 


0(6*)= P(- cos#) = P(w) 


(B20) 


(The choice of w = cos 0 would also be acceptable and indeed is the most popular one. 
However, employing Eq.(B20) allows the solution to be regular at 0 = n and to contain a 
singularity at 0 = 0 ; this property will be essential when matching source boundary 
conditions to the problem.) The solution to Eq.(B19) will be taken to be given by the 
infinite series in w, i.e., 


P{w)=Y J a j w i (B21) 

j = o 

which gives, upon substitution into Eq.(B19), 


Z(/ + 2 \j + tyj+2 wJ ~ Z./(/ + tyjW J + ^Ej a j wJ = 0 (B22) 

j= o j= o j= 0 

Since this relation must hold for every power of w, the coefficient separately must satisfy 

(J + 2 XJ+ l ) a . /+ 2 - ./(/ + 1 ) a j + to J = 0 (B23 ) 

Solving this equation for a j+2 yields 


X/ +1 )-^ 

(i + 2X7 + 1)° 


(B24) 
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If the series given by Eq.(B21) does not terminate at some value for j , it is seen from 
Eq.(B24) that 


lim^^- = — — (B25) 

7-*° a, ; + 2 

J u 

At the values of w = ±1 , the series will diverge. This is not an acceptable situation for a 
physically realizable situation. Thus, the series must terminate at some value of j , e.g., 
j = 1. From Eq.(B24), with a l+2 = 0 (i.e., the series tenninates), one must have 

A = 1(1+ 1) (B26) 

where / is a positive integer. This is the fundamental source of the combination /(/ + 1) 
throughout the equations of propagation involving the angular coordinate 9. Therefore, 
Eq.(B19) becomes 


_d_ 

dw 


(i-+) 


dP } (vv) 
dw 


+ l(l+l)P,(w) = 0 


(B27) 


in which P,(w) is now also a function of /. This differential equation defines [1] the 
Legendre functions P,(w). 

To find explicit expressions for the functions P,(w), consider another change of 
variables, i.e., 


o = 



(B28) 


Equation (B27) then simplifies to 


(l-w i y^- + 2lwu = 0 

V J dw 


Differentiating this relation with respect to w j + I times gives 

(! - w2 )^J- & ~ 2l + 2 ^ W ^, + ( 2l + = 0 


(B29) 


(B30) 


where 


d'o 

t = r 

dw' 



At the specific value of j = /, Eq.(B30) yields 


(B31) 
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(B32) 


(i -w 1 y^-^--2w-^-+ /(/ +i y = o 

v >dw 2 dw v ’ 

which is Eq.(B27) for t = P ; (w). At j = 1, Eq.(B31) differs only be a constant factor. 
Thus, at the outset, setting E’ ; (l) = l for all l finally gives 

PXw) = ^r(w 2 -lj (B33) 

’ 2 l l\dw ,y ' 

The first three Legendre functions are thus given by 

nW=l (B34) 

P x (w)=w = -cosO (B35) 

P 2 (w)= ±(3w 2 - 1)= ^(3cos 2 0- l) (B36) 


The analysis of Eq.(B7) can proceed in one of two ways. First, it must be noted 
that like that of the function P t (w ), the radial function now also becomes a function of / 
upon use of Eq.(B26), i.e., i?(r) = R,(r\ In the first approach to Eq.(B7), one can simply 
expand the second derivative with respect to r and obtain 


d~R, 2 dR, 

— y + — - 

dr r dr 


+ 


1(1+ ij 


2 

r j 


R,= 0 


(B37) 


Solutions to this equation are linear combinations of spherical Bessel Functions [2,3] of 
the first kind j,(kr) and of the second kind n l (kr ), i.e., 

R, (r) = Ajj (kr) + Bn , (kr) (B 3 8 ) 

where the constants A and B are detennined by the boundary conditions of the problem. 
The analysis of these functions can be carried out in an analogous manner to that of 
Eq.(B19); see, e.g., [3] for details. 

Hence, for the propagation problem addressed in this work, the solution of 
Eq.(Bl) is given by Eqs.(B3), (BIO), (B 1 7), and (B20), viz, 

u e (r, e) = R, (r)P, (~ cos 9) (B39) 

or some linear combination of this composite function, adhering to some prevailing 
boundary and source conditions. 
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A second approach using Eq.(B7) in terrestrial propagation scenarios involves the 
use of a transformation and an approximation. First, the transformation simply involves 
the prescription 


Ri(r)=rR,(r) 


(B40) 


allowing Eq.(B7) to become 




dr 2 


k- 


r 2 J 


R,= 0 


(B41) 


The general solution of Eq.(B39) now becomes in tenns of the radial function R,(r) 


U e (r, 0)=C P/ (- cos 0) 


(B42) 


Elsing the solution of Eq.(B41) is, remembering Eqs.(B38) and (B40), 

R,(r)= r\_Aj,{kr)+ Bn,(kr )] (B43) 


from which, using Eq.(B42), leads back to Eq.(B39). The point here is that the form of 
Eq.(B41) is amenable to approximations which give solutions that differ from Eq.(B43). 
This fonn of the radial equation is used in this work because of this flexibility. An 
example of such an approximate solution will now be given. 

Application of Eq.(B41) to a terrestrial propagation problem above the surface of 
a planet of radius a is as follows. Letting E, be the distance above the surface, one has 
r = a + where, in general, g « a . In this instance, one can approximately write 


± = (a + £f = 


a 1 

( £] 

-2 

« a 1 

( 

1-2— 


l aJ 


f aJ 


Using this in Eq.(B41) and changing variables from r to Eg yields 

d-R(S) 


tli- 


■ + 


A '(/+!), 2/(/+ 1) 2 


cr ) 




Changing variables once again using 




^ 2/(/ + l) V/3 

V 7 


4 


(B44) 


(B45) 


(B46) 


finally gives 
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(B47) 


d%(£) , 

df 2 



/ ( /+1 ) 


v 


a 2 A2/(/+l). 


2/3 


+ C 


^,(0=o 


Letting 


f 



l(l+l) 


V 


2 

a A 


\ 2/3 


2i(j + i)J 


Eq.(B47) can be written as 




[f,-£]*,(0 = O 


(B48) 


(B49) 


the solution of which are Airy functions [4], This form of the radial equation was 
extensively employed by Fock [5] in early investigations of radio wave propagation 
above the Earth’s surface; see also [6]. The specific values of t, are determined by the 
prevailing boundary conditions. 

Finally, the expressions for the boundary conditions for the tangential electric and 
magnetic field discussed in Appendix I, (Eqs.(A48) and (A50)) can be expanded upon. 
Equation (A48) for the continuity of the tangential electric field across a boundary at 
r = a separating regions 1 and 2 becomes, using Eq.(B42), 


AAl 

3* 


dr 


(B50) 


Similarly, Eq.(A50) for the continuity of the tangential magnetic field across the 
boundary becomes 




= n l \K: 


(B51) 
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